Crustal Entrainment and Pulsar Glitches 
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Large pulsar frequency glitches are generally interpreted as sudden transfers of angular momentum 
between the neutron superfiuid permeating the inner crust and the rest of the star. Despite the 
absence of viscous drag, the neutron superfiuid is strongly coupled to the crust due to non-dissipative 
entrainment effects. These effects are shown to severely limit the maximum amount of angular 
momentum that can possibly be transferred during glitches. In particular, it is found that the 
glitches observed in the Vela pulsar require an additional reservoir of angular momentum. 

Keywords: neutron star, pulsar glitches, superfluidity, entrainment 



Introduction. Since their fortuitous discovery by Joce- 
lyn Bell and Anthony Hewish in 1967, more than 2000 
pulsars have been found [1]. Their identification as neu- 
tron stars the compact residues of type II supernova 
explosions predicted by Baade and Zwicky in 1933 fsj, 
was definitively established the next year after the dis- 
coveries of pulsars in the Crab and Vela supernova rem- 
nants. Pulsars are among the most accurate clocks in 
the universe with periods ranging from 1.4 milliseconds 
up to several seconds. The delays associated with the 
spin-down of the star are at most of a few of milliseconds 
per year. 

Nevertheless, irregularities have been detected in long- 
term pulsar timing observations f^]. In particular, some 
pulsars exhibit sudden increases in their rotational fre- 
quency ri. These "glitches" , whose amplitude varies from 
Afl/n ~ 10~^ up to ~ 10~^ are generally followed by a 
relaxation over days to years and are sometimes accom- 
panied by a sudden change of the spin-down rate from 
An/n\ - 10"*^ up to - 10^2 (ggg g Section 12.4 in 

!)• 

Soon after the first observations of glitches in the Vela 
and Crab pulsars, several scenarios were advanced @]. In 
particular, glitches were thought to be the manifestations 
of starquakes, but this could not explain the frequent oc- 
currence of Vela pulsar glitches 0] ■ A corequake model of 
Vela pulsar glitches was proposed [l2|, but the existence 
of a solid core later appeared to be highly speculative 
(see e.g. 0). The long relaxation times following glitches 
provided strong evidence for the presence of superfluids 
in neutron star interiors and hinted at its possible role 
in the glitch mechanism itself P, ^ . Neutron-star super- 
fluidity had been predicted and studied even before the 
discovery of pulsars [Hj. Anderson and Itoh devel- 
oped the fruitful idea that Vela like glitches are related to 
the dynamics of the neutron superfluid permeating the 
inner crust of neutron stars [ist . 

Vortex-mediated glitches. The neutron superfluid is 
weakly coupled to the crust by mutual friction forces and 
thus follows its spin-down via a radial motion of quan- 
tized vortices away from the rotation axis unless vortices 
are pinned to the crust. In this case, the superfluid can 



rotate more rapidly than the crust. The lag between the 
superfluid and the crust induces a Magnus force acting 
on the vortices thereby producing a crustal stress. When 
the lag exceeds a critical threshold, the vortices are sud- 
denly unpinned. As a result, the superfluid spins down 
and, by the conservation of the total angular momentum, 
the crust spins up leading to a glitch. This scenario found 
some support from laboratory experiments in superfluid 



helium [IJ, [15| . The good fit to the glitch data triggered 



further developments to explain the postglitch relaxation 
by the motion of vortices lit- isf. 

In the meantime, it was argued that the core (sup- 
posed to contain superfluid neutrons and type I super- 
conducting protons) is unlikely to play any role in glitch 
events |l^ (see also (2^). Due to non-dissipative entrain- 
ment effects similar to those arising in superfluid '^He-'''He 
mixtures, neutron superfluid vortices carry a fractional 
magnetic quantum flux. Electron scattering off the mag- 
netic field of the vortices leads to a (dissipative) mutual 
friction force acting on the superfluid. As a result, the 
core superfluid is strongly coupled to the crust and to the 
charged particles, thus following the long-term spin-down 
of the star caused by electromagnetic radiation. 

The confidence in the vortex-mediated glitch inter- 
pretation led to a new constraint on the structure of 
neutron stars hence also on the equation of state of 
dense matter 12111 . The latest models like the "snow- 
plow" model [22| can reproduce various observations of 
pulsar glitches. However, many fundamental aspects of 
these models remain poorly understood. For instance, 
the strength of vortex pinning, which is one of the cru- 
cial microscopic inputs, has been a controversial issue 
over the past years (see e.g. Section 8.3.5 of Q). The 
mechanism that triggers the unpinning of vortices like 
superfluid instabilities is also a matter of debate. 

More importantly, these models ignore the non- 
dissipative entrainment effects in neutron-star crusts 
which have been shown to be very strong [H-iil. In 
this Letter, the impact of crustal entrainment on pulsar 
glitches is studied. Using the latest pulsar glitch data Q , 
it is shown that the neutron superfluid in the crust does 
not carry enough angular momentum to explain the Vela 
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pulsar glitches. 

Entrainment in neutron-star crusts. It has been re- 
alized only recently that entrainment arises not only in 
the core of a neutron star but also in the crust because 
unbound neutrons can be elastically scattered by the 
crustal lattice for specific wave vectors, as determined 
by Bragg's law 24-2^. A neutron that is Bragg re- 
flected cannot propagate and is therefore entrained by 
the crust. Unlike viscous drag, entrainment is non- 
dissipative. Even if a neutron is not Bragg reflected, its 
motion will still be affected by the crustal lattice. Neu- 
tron diffraction experiments are routinely performed to 
study crystal structures. The specificity of neutron-star 
crusts is that neutrons form a highly degenerate quantum 
liquid. Due to the Pauli exclusion principle, neutrons 
have different wave vectors and are therefore diffracted 
differently. Entrainment can be characterized by the den- 
sity of conduction neutrons, i.e. neutrons that are 
effectively "free" , or equivalently by an effective neutron 
mass m* = TO„n^/n^ where m„ is the bare neutron mass 
and the density of unbound neutrons. Neutron con- 
duction has been systematically studied in all regions 
of the inner crust using a state-of-the art crust model 
based on the band theory of solids [i^ . Entrainment has 
been found to be very strong, especially in the intermedi- 
ate region of the inner crust at average baryon densities 
h ^ 0.02 — 0.03 fm^"^ where ^ or equivalently 

TO* J> TO„. 

Pulsar glitch constraint. Large pulsar glitches are usu- 
ally interpreted as sudden transfers of angular momen- 
tum between the neutron superfluid in the crust and the 
rest of the star 17|. This model predicts that the ra- 



with 



tio of their respective moments of inertia must obey the 
constraint 
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(1) 



where Ag is the glitch activity parameter defined by the 
sum over glitches occurring during a time t 



A„ = 



(2) 



while Q, is the average spin-down rate. Both Ag and 
can be measured from pulsar-timing observations. Since 
Is <C Ic, Ic can be replaced by the moment of inertia 
I = Is + Ic of the entire star. Approximating Is by the 
moment of inertia /crust of the crust, a constraint on the 



mass and radius of the Vela pulsar was derived in [21 1. 
This approximation treats all unbound neutrons as con- 
ducting {n'^ = nl), an assumption which turns out to 
be unrealistic 24, 2^ 2^. Due to entrainment, the an- 
gular momentum Jg of the superfluid depends not only 
on the angular velocity fls of the superfluid, but also on 
the observed angular velocity 51 of the star and can be 
expressed as [lO] 

Js = Iss^s + (Is - Iss)^ , (3) 



Is = J mnul^g^d^r , Iss ^ ym^n^p^d^- 



(4) 



where g is the cylindrical radius. The constraint ([T} thus 
becomes 1301 
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This inequality is much more stringent than ([T]) because 

/ss > Is- 

Results. The ratio appearing in the left hand side of 
Eq. ([5]) can be decomposed as 
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In the thin crust approximation 31], Iss is given by 
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where M and R are the neutron-star mass and radius, p 
is the average mass density, Pcorc is the pressure at the 
crust-core transition and -Pdrip(^ /core) is the neutron- 
drip pressure. Corresponding expressions for Is and /crust 
are obtained by replacing nl^m^ in Eq. ([7]) by n\jnn and 
p respectively. Note that /ss //crust and /s //crust depend 
only on the crust properties, and can be written as 
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where n = pjmn is the average baryon density. Integrat- 
ing Eqs. (HI) and ^ with the trapezoidal rule using the re- 
sults of ^29j| summarized in Table HI we find Iss — 4.6/crust 
and Is ~ 0.89/crust leading to (/s)^//ss ^ 0.17/crust- The 
ratio /crust // depends on the global structure of neutron 
stars. We have made use of Eq. (47) of 3^. This formula 
was obtained by solving the equations of General Rela- 
tivity using a set of realistic dense-matter equations of 
state (EoS). Results for (/s)^/(//ss) are shown in Fig[T] 
Note that microscopic calculations based on chiral effec- 
tive field theory (33| (and more generally any realistic EoS 
as illustrated in Fig. [3|) as well as observations of X-ray 
binaries indicate that neutron stars with M = Mq 
have a radius i? ^ 13 km. 

Because it was the first observed pulsar to exhibit 
glitches. Vela has become the testing ground for glitch 
theories. Since 1969, 17 glitches have been detected i4j]. 
As shown in Fig.O the cumulated glitch amplitudes given 
by AQi/Q ~ tAg (with an appropriate choice of time 
origin) increases almost linearly with the time t. A lin- 



ear fit yields Ag ~ 2.25 x 10" 



With the angular 
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P (MeV fm"^) 


ni/n (%) 


n-n/ni (%) 


0.0004575 


15.0 


82.6 


0.0009886 


61.1 


27.3 


0.006097 


82.6 


17.5 


0.01507 


86.0 


15.5 


0.03820 


87.9 


7.37 


0.06824 


89.1 


7.33 


0.1068 


86.6 


10.6 


0.1561 


89.1 


30.0 


0.2183 


89.2 


45.9 


0.2930 


89.4 


64.6 


0.3678 


100 


64.8 



TABLE I. Entrainment parameters in the inner crust of cold 
non-accreting neutron stars as obtained in [2^. n is the av- 
erage baryon density, is the density of unbound neutrons, 
and nJi is the density of conduction neutrons. The pressure 
P was calculated using the formulae in Appendix B of [sj- 
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FIG. 1. (Color online) (Ib)'^ / (I Ibb) for different neutron-star 
radii R and masses M from IMq (upper curve) to 2Mq (lower 
curve). The shaded area is excluded if Vela pulsar glitches 
originate in the inner crust. 
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FIG. 2. Cumulated glitch amplitudes as a function of the 
modified Julian date for the Vela pulsar from Q (square 
symbols) and linear fit (solid line). 
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FIG. 3. (Color online) Neutron-star mass-radius diagram for 
three different unified equations of state. The shaded area is 
excluded by Vela pulsar glitch data, assuming that only the 
neutron superfluid in the crust is involved. The red (blue) 
area is the constraint obtained with (without) taking into 
account crustal entrainment. The sensitivity of these con- 
straints with respect to the crust-core transition pressure Pcorc 
is indicated by the dashed line and the dotted line. 



frequency Q — 11.1946499395 Hz and average spin-down 
rate fl = -1.5666 x 10"" s'^ we find G ~ 1.6%. A 
similar estimate has been obtained from a statistical anal- 



ysis of Vela like pulsars [35|. As illustrated in Fig.[Tl com- 
bining the glitch data with Eq. ([S]) leads to a constraint 
on the global neutron-star structure. This constraint, 
which can be approximately written as i? > 8.51-I-5.23M 
or cquivalently M < 0.190i?- 1.61 with M in Mg and R 
in km, is plotted in Fig. |31 together with three represen- 
tative unified EoS spanning different degrees of stiffness 
of dense neutron matter, from the softest (BSkl9) to the 
stiffest (BSk21), as obtained from microscopic calcula- 
tions [36|. The sensitivity of the glitch constraint with 
respect to the corresponding crust-core transition pres- 
sure is also shown. This analysis implies that Vela (and 
more generally pulsars with Vela like glitches) should be 
less massive than our Sun {AI < 0.6Mq for the soft- 
est EoS). Such a low mass neutron star is unlikely to be 
formed in a type II supernova explosion f38|. However, 
the association of the Vela pulsar with the eponymous 
supernova remnant is well established. 

Discussion. We are thus led to conclude that the 
neutron superfluid in neutron-star crusts does not carry 
enough angular momentum to explain large pulsar 
glitches like those observed in Vela, unless crustal en- 
trainment and crust-core coupling are much weaker than 
considered here. A similar conclusion has been reached 
in d^. 

The estimates of m* obtained in [2^ agree closely with 



previous calculations [24| using a different model thus 
suggesting that strong crustal entrainment is generic. 
Moreover, was found to be weakly dependent on 
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the crystal structure 2J]. The existence of nuclear 
"pasta" phases near the crust bottom (see e.g. Sec- 
tion 3.3 of 0]), which we have ignored, might enhance 
the neutron conduction owing to the low dimensional- 
ity of such configurations. However, it has been ar- 
gued that these pastas (if any) could only exist in a 
small region of the crust, at baryon densities above 
~ 0.06 fm~'^, if the lowest-frequency quasiperiodic os- 
cillation observed in giant flares from soft gamma-ray 
repeaters (SGRs) is to be interpreted as the fundamen- 
tal torsional crustal mode Setting = nl^ for 
n > 0.06 fm~^, the impact of pastas is found to be 
small since the ratio Iss/Icmst is reduced from 4.6 to 4.3 
whereas (/s)^/(^rust-^ss) is raised from 0.17 to 0.19. In 
reality, is never equal to in any region of the crust, 
even in the presence of pastas On the other hand, 
the spin-orbit coupling (which was neglected in 12911 1 . 
would increase the number of entrained neutrons [2J] 
and could be more important than pastas since it op- 
erates at all densities. We also anticipate that quantum 
and thermal fluctuations of ions about their equilibrium 
positions, crystal defects, impurities and more generally 
any kind of disorder would presumably reduce (but not 
cancel) entrainment effects. Further work is needed to 
confirm these espectations. 

On the other hand, the strong coupling of the core to 
the crust that we have considered here relies on the as- 
sumption of type I superconductivity [3| ■ The observed 
rapid cooling of the neutron star in Cassiopeia A has 
recently provided strong evidence for core neutron su- 
perfiuidity and proton superconductivity, but not on its 
type [m ]^ . If the superconductor is of type II, the cou- 
pling time could be much longer [43[. On the other hand, 
type II superconductivity has not only been argued to be 
incompatible with observations of long-period precession 
in pulsar s |44| but has also been questioned on theoretical 
grounds [45j. In fact, the superconductor might neither 
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be of type I nor of type II 46[ . In addition, neutron-star 



cores might contain other particle species with various 
superfluid and superconducting phases. 

Removing the discrepancy between glitch models and 
observations thus requires a closer examination of crustal 
entrainment and crust-core coupling. The regularity of 
glitches illustrated in Fig. [5] and the fact that G ^ 2% 
suggest the existence of a reservoir of angular momentum 
in a limited region of the star, possibly in the outermost 
part of the core just below the crust (e.g. [13, li^). This 
warrants further studies. 

This work also shed light on the importance of crustal 
entrainment, which has been generally overlooked even 
though they may have implications for other astro- 
physical phenomena like quasiperiodic oscillations in 
SGRs Hi. 
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